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Abstract We have studied the structure of hot accretion 
flow bathed in a general large-scale magnetic field. We have 



considered magnetic parameters 



where z are the Alfven sound speeds in three direction 
of cylindrical coordinate (r, ip, z). The dominant mecha- 
nism of energy dissipation is assumed to be the magnetic 
diffusivity due to turbulence and viscosity in the accretion 
flow. Also, we adopt a more realistic model for kinematic 
viscosity [y = ac s H), with both c s and H as a function 
of magnetic field. As a result in our model, the kinematic 
viscosity and magnetic diffusivity (77 = r]QC s H) are not con- 
stant. In order to solve the integrated equations that govern 
the behavior of the accretion flow, a self-similar method is 
used. It is found that the existence of magnetic resistivity 
will increase the radial infall velocity as well as sound speed 
and vertical thickness of the disk. However the rotational 
velocity of the disk decreases by the increase of magnetic 
resistivity. Moreover, we study the effect of three compo- 
nents of global magnetic field on the structure of the disk. 
We found out that the radial velocity and sound speed are 
Sub-Keplerian for all values of magnetic field parameters, 
but the rotational velocity can be Super-Keplerian by the 
increase of toroidal magnetic field. Also, Our numerical 
results show that all components of magnetic field can be 
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important and have a considerable effect on velocities and 
vertical thickness of the disk. 

Keywords accretion, accretion disks, magnetic field, mag- 
netohydrodynamics:MHD 



1 Introduction 

The study of advection accretion flows around low lumi- 
nosity black hole candidates and neutron stars is currently a 
very active field of research, both theoretically and observa- 
tionally (for a review, see Narayan, Mahadevan & Quataert 
1998). Observational evidence for the existence of low- 
luminosity black holes at the centre of galaxies and in AGNs 
(Ho 1999) makes it necessary to revise theoretical models of 
the accretion discs. Thereby, the development of the subject 
of advection-dominated accretion flows (ADAFs) in recent 
years has led to global solutions of advection accretion discs 
around accreting black hole systems and neutron stars. In 
this case, viscously generated internal energy is not radiated 
away efficiently as the gas falls into the potential well of the 
central mass (as in standard thin disc models; Shakura & 
Sunyaev 1973). Instead, it is retained within the accreting 
gas and advected radially inward (Narayan & Yi 1994, here- 
after NY 1994) and might eventually be lost into the central 
object or, in contrast, a considerable portion of it might give 
rise to wind on to black holes and neutron stars (Blandford 
& Begelman 1999). By definition, ADAFs have very low ra- 
diative efficiency and as a consequence they can be consid- 
erably hotter than the gas flow in standard thin disc models 
(Narayan & Yi 1995a,b); therefore, they are ultra-dim for 
their accretion rates (Phinney 1981; Rees et al. 1982). 

A notable problem arises when the accretion disc is 
threaded by a magnetic field. In the ADAF models, the tem- 
perature of the accretion disc is so high that the accreting 
materials are ionized. The magnetic field therefore plays an 
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important role in the dynamics of accretion flows. Some au- 
thors have tried to solve the magnetohydrodynamics (MHD) 
equations of magnetized ADAFs analytically. The effect of 
toroidal magnetic field on the disc were studied for exam- 
ple by Akizuki & Fukue 2006, Abbassi et al 2008, 2010 and 
Faghei 2012. Also Zhand & Dai 2008 have already consid- 
ered global magnetic field on the disc. Ghanbari et al 2007, 
2009 present the effect of viscous-resistive in ADAFs bathed 
in a dipole magnetic field. Resistive diffusion of magnetic 
field would be important in some systems, such as the pro- 
tostellar discs (Stone et al. 2000; Fleming & Stone 2003), 
discs in dwarf nova systems (Gammie & Menou 1998), the 
discs around black holes (Kudoh & Kaburaki 1996), and 
Galactic center (Melia & Kowalenkov 2001; Kaburaki et al. 
2010). In ADAF models, energy dissipation in the accre- 
tion flow can be assumed to be a result of turbulent viscosity 
and electrical resistivity. The Magnetic energy density of the 
flow must be dissipated by ohmic heating with a rate compa- 
rable to that of the viscous dissipation (Bisnovatyi-Kogan & 
Ruzmaikan 1974, Bisnovatyi-Kogan, Lovelace 1997). We 
argue that the ohmic and viscous dissipation must occure as 
a result of plasma instabilities. Bisnovatyi-Kogan, Lovelace 
1997 point out that ohmic dissipation of magnetic energy 
density has an important role in heating process if accretion 
flows with condition of equipartition, e mag ~ £kin> where 
£mag&£kin are the magnetic and kinetic energy density, re- 
spectively. In this regard, note that although Narayan & Yi 
(1995) assume an equipartition magnetic field, they do not 
consider the ohmic dissipation. 

Under some conditions, it is important that we consider 
the effect of resistivity on accretion flows. Kuwabara et al. 
(2000) showed the results of global magnetohydrodynamic 
(MHD) simulations of an accretion flow initially threaded by 
large-scale poloidal magnetic fields, including the effects of 
magnetic turbulent diffusivity. They found the importance 
of the strength of magnetic diffusivity when they studied it 
in magnetically driven mass accretion. They pointed out that 
the mass outflow depends on the strength of magnetic diffu- 
sivity, so that for a highly diffusive disc, no outflow takes 
place. 

In this paper, we extend the work of Akizuki & Fukue 
2006 and ghanbari et al. 2007, Ghanbari et al 2009 and 
Faghei 2012 by considering large-scale magnetic field with 
three components in cylindrical coordinates (r, ip, z) in a 
viscose-resistive accretion disc, and investigate the role of 
non-constant magnetic diffusivity in the system. Ghan- 
bari, Salehi & Abbassi (2007) have presented a set of self- 
similar solutions for two-dimensional (2D) viscous-resistive 
ADAFs in the presence of a dipolar magnetic field of the 
central accretor. They have shown that the presence of a 
magnetic field and its associated resistivity can considerably 
change the picture with regard to accretion flows. The paper 
is organized as follows. In section 2 we present the basic 



magnetohydrodynamics equations, which include the three 
components of magnetic field and magnetic resistivity. Self- 
similar equations are investigated in section 3 and the sum- 
mery will come up in section 4. 



2 Basic Equations 

Since we are interested in analyzing the structure of a mag- 
netized ADAF bathed in a global magnetic field, we con- 
sider a magnetic field in the disc with three components, 
B r , and B z . we suppose that the gaseous disc rotating 
around a compact object of mass M*. Thus, for a steady axi- 
symmetric accretion flow, i.e., d/dt = d/dip — , we can 
write the standard equations in the cylindrical coordinates 
(r, cp, z). We vertically integrated the flow equations and, 
all the physical variables become only function of the ra- 
dial distance r. Moreover, we neglect the relativistic effects 
and Newtonian gravity in radial direction is considered. For 
conservation of energy, it is assumed the energy which is 
generated due to viscosity and resistivity dissipation are bal- 
anced by the radiation and advection cooling. Under these 
assumption, we can rewrite the MHD equation as (Zhang & 
Dai 2008). So the equation of continuity gives 



ld_ 

r dr 



(rl2v r ) = 2pH, 



(1) 



where E is the surface density at the cylindrical radius r, 
which is define as £ = 2pH, v r the radial infall velocity, 
p the mass loss rate per unit volume, H would be the disc 
half-thickness and p is the density of the disc. 
The equation of motion in the radial direction is 



dv r id GM* 1 d 



1 d 



dr 



r r* 
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where v v is the rotational velocity, c s the isothermal sound 
speed, which is define as c 2 = p gas /p, Pgas being the gas 
pressure. Here, c r , c v and c z are Alfven sound speeds in 
three direction of cylindrical coordinate and are defined as 



R 2 



(3) 



The equation of angular momentum transfer can be writ- 
ten as 



v r d 
r dr 



1 d , , „df2, c r d , , — 



(4) 
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where f2(= v^/r) is the angular velocity. Also, we assumed 
that only the (ys-component of viscose stress tensor is impor- 
tant which is t rtp = firdfl/dr , where p,(= up) is the vis- 
cosity and v is the kinematic coefficient of viscosity. 

The hydrostatic balance in vertical direction is integrated 

to 



n 2 K H 2 - -L Cr -?-(VZc z )H = c 
VS dr 



ci), 



(5) 



Now we can write the energy equation considering cool- 
ing and heating processes in an ADAF. Therefore the energy 
equation will be 



del 



7 



1 dr 



v r c Q 



,dp 
dr 



Q+-Q- 



(6) 



where 7 is the radio of specific heats. In the right hand side 
of energy equation, Q + = Q V i S + Qb is the dissipation rate 
by viscosity Q v i s and resistivity Qb, and also Q_ = Q ra d 
represent the energy lose through radiative cooling. For the 
right hand side of energy equation we can write 



Qadv — Q+ — Q- — fQ-i 



(7) 



Here, Q a dv represents the advection transport of energy and 
is defined as the difference between the magneto-viscose 
heating rate and radiative cooling rate. We employ the ad- 
vection parameter / = 1 — to measure the hight de- 
gree to which accretion flow is advection-dominated. When 
/ ~ 1 the radiation can be neglected and accretion flow is 
advection dominated while in the case of small /, disc is in 
the radiation dominated case. Now the viscous and resistive 
heating rate are expressed as 



Q 



vpr (— ) 
dr 
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(8) 
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where J = V x B is the current density. We assume both 
the kinematic viscosity coefficient v and the magnetic diffu- 
sivity f] have the same unites and to be due to turbulence in 
the accretion flow. Bisnovatyi-Kogan & Lovelace 1997 have 
also shown that the flows have an equipartition magnetic 
field with the results that dissipation of magnetic energy at 
a rate comparable to that turbulence must occur by ohmic 
heating. They argue that this heating occurs as a results 
of plasma instabilities. Then, it is physically reasonable 
to express 77 such as v via the a-prescription of Shakura & 
Sunyaev (1973) as follow (Bisnovatyi-Kogan & Ruzmaikin 
1976) 



ac s H 



(10) 



77 = r/ c s H 



(ID 



where, a and 770 are the standard viscous parameter and the 
magnetic diffusivity parameter respectively. Also these pa- 
rameters are assumed to be positive constant and less than 
unity (a, 770 < 1) (Compbell 1999; Kuwabara et al. 2000; 
King et al. 2007). 

Finally since we consider a global magnetic field, the 
three components of induction equation can be written as: 



B r = 



B v = — [{v v B r - v r B v ) - -irB v )\ 



1 d r . dB z , 

— ^[r(v r B z + V ^)] 



dr 



(12) 



(13) 
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where B r . v . z are the field scaping/creating rate due to mag- 
netic instability or dynamo effect. Now we have a set of 
MHD equations that describe the structure of magnetized 
ADAFs. The solutions to these equations are strongly corre- 
lated to viscosity, resistivity, magnetic field strength P r ,tp,z 
and the degree of advection / . We seek a self-similar so- 
lution for the above equations. In the next section we will 
present self-similar solutions to these equations. 



3 Self-Similar Solutions 

In order to have a better understanding of the physical pro- 
cesses taking place in our discs, we seek self-similar solu- 
tions of the above equations. The self-similar method has a 
wide range of applications for the full set of MHD equations 
although it is unable to describe the global behavior of ac- 
cretion flows since no boundary conditions have been taken 
into account. However, as long as we are not interested in 
the behavior of the flow near the boundaries, these solutions 
are still valid. In the self- similar model the velocities are 
assumed to be expressed as follow 



v r (r) 



-C\OL 



GM* r 
ro r 



!GM*.r . 

c 2\l ( — )" 

ro r a 



^(r) = c 3 (^)(^)- 1 



^0 r 



r>2 



GAL 



(15) 
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where constants c\,c<i and C3 are determined later from the 
main MHD equations. Also ro are exploited to write the 
equation in non-dimensional form and the constants f3 r ,v,z 
measure the radio of the magnetic pressure in three direction 
to the gas pressure, i.e., @ r ^ z = p mag ,r, v ,z/Pgas- 
Assuming the surface density S to be in the form of 



E(r) = £„(-)' 
ro 

where Eq and s are constant, we obtained,e.g. 



p = po(-r 



B 



.2 — BrO.tpO.zoi — ) 2 

ro 



(20) 



(21) 



where p and B r o tV o,zO are constants. The half-thickness of 
the disc still satisfies the relation H oc r and we can write 



H(r) = c 4 r (-) 

7-0 



(22) 



By substituting the above self-similar solutions in the conti- 
nuity, momentum, angular momentum, hydrostatic balance 
and energy equation of the disc, we obtain the following sys- 
tem of dimensionless equations to be solved for ci,C2, C3 and 
c 4 : 



lsdCiEo /GM, 
is + -) 



2 ' 2ric 4 



ro 



(23) 



? = c§-l-[(*-l) + (*-l)&+(*+l))8 v ]ca (24) 



- 22 2 
-a c - 



-acic 2 = -2^ s + l) ac 2\/c3C4 + (s + 1)c 3 a/5^^ (25) 



c 4 = ^ - l) 2 /3r/Vs + 4(1 + J r + /3 v )c 3 



+ (a - l)V^Cs] (26) 



1 9 
(—j - ( s - l))aciv^ = ^/ac2C 4 

+ ^/%c 3 C4 + (s - 2) 2 /?,] (27) 

If we solve the self-similar structure of the magnetic field 
escaping rate, we obtain 



Bor = 0, 



(28) 
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47r/3 y c 3 S]o 
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(19) Bo, = ( : 



47r/3 z c 3 E 
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1 GM* 

> 5/2 



It is evident from the above equations that for s = —1/2 
there is no mass loss while it is present for the case where 
s > —1/2. On the other hands, the toroidal component of 
escape and creation of magnetic fields balance one another 
for s — 3 and also for z-component of field creating / escap- 
ing rate ("eqf30b they balance for s = 1. Although outflow is 
one of the most important processes in accretion theory, (see 
Narayan & Yi 1995; Blandford & Begelman 1999; Stone & 
Pringle & Begelman 1999 and also some recent work like 
Xie & Yua 2008; Ohsuga & Mineshige 2011), however in 
our model we choose a self-similar solution in which p = 0, 
(s = —1/2) thus ignoring the effect of wind and outflow on 
the structure of the disc. 

Without magnetic resistivity and magnetic field, 77 = 
Pr = Ptp = Pz = 0, the equations and their similarity solu- 
tions are reduced to the standard ADAFs solution (Narayan 

6 Yi 1994). Also without resistivity they are reduced to 
Zhang & Dai 2008. 

Now, we can do a parametric study considering our input 
parameters. The parameters of our model are the standard 
viscose parameter a, the magnetic diffusivity parameter 770, 
the advection parameter /, the radio of the specific heats 

7 and the degree of magnetic pressure to the gas pressure 
in three dimensions of cylindrical coordinate, (3 r , f3 v and 
p z . Figure Q] shows how the coefficients ac±(= —v r /i)K, 
ci{= v v /v K ), c s/vk) and c 4 (= H/r) depend on 
the magnetic parameter in radial direction f3 r for several val- 
ues of magnetic diffusivity parameter r/ , i. e. t/q = (dot- 
ted line) , 770 = 0.05 (dashed line) and 770 = 0.1 (solid line), 
corresponding to a = 0.1, 7 = 4/3, /3 V = (3 Z = 1.0 and 
/ = 1.0 (fully advection). radial velocity is shown in the 
upper-left panel. In ADAFs the radial velocity is generally 
less than free fall velocity on a point mass, but it becomes 
larger if the magnetic parameter f3 r is increased. Also, we 
can see that the radial flows of the accretion materials be- 
come larger by the increase of resistivity parameter 770. The 
upper right panel of figure[T]shows the radio of the rotational 
velocity to the Keplerian one. As it can seen, the rotational 
velocity slightly shifts up when the magnetic parameter f3 r 
increases. Although the increase of resistivity parameter 7/0 
will decrease the rotational velocity. Moreover, the down 
left and right panels of figure[T]display the radio of the sound 
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Fig. 1 Numerical coefficient c^s as a function of magnetic parameter j3 r for several values of 770. The dotted, dashed and solid lines 
correspond to rjo = 0.0, 0.05 and 0.1 respectively. Parameters are set as s = — 0.5(no wind), a = 0.1, = /3 Z = 1 and f = 1. 



speed to the Keplerian velocity and the vertical thickness of 
the disc respectively. As the magnetic parameter j3 r become 
large, the isothermal sound speed of the flow decreases but 
vertical thickness of the disc has a different behavior be- 
tween /3 r < 1 and f3 r > 1. When /3 r is below unity, the 
vertical thickness is reducing. On the other hand, for /3 r > 1 
the thickness is increasing. It means that there is a mini- 
mum near j3 r ~ 1. Also, by adding the magnetic diffusivity 
parameter, the sound speed and vertical thickness will in- 
crease. 

In figure |2]behavior of the coefficients aci, C2, C3 and C4 
versus toroidal magnetic field B v (j3 v ) are shown for dif- 
ferent values of r)Q. From the upper-left panels of figure [2] 
we can see that for f3 v = — 10, the radial infall velocity 
is sub-Keplerian, and that becomes larger by increasing of 
resistivity parameter t/q. In addition, when the toroidal mag- 
netic field becomes stronger, the radial velocity of accretion 
materials increases. The upper-right panel in figure ^dis- 
play the rotational velocity of accretion disc. We see that for 
the given advection parameter /(= 1.0) and f3 r = /3 Z = 1.0, 
the rotational velocity increase as the toroidal magnetic field 
become stronger. The considerable matter about this panel 
is that, when the toroidal magnetic field is very strong, the 
rotational velocity will be Super-Keplerian (i.e., /3„ > 7). 
Also as the magnetic diffusivity parameter rj increases, the 



rotational velocity of the disc decrease. The raise of viscose 
torque by adding 770 and /3 V parameters generate a large neg- 
ative torque in angular momentum equation and cause the 
angular velocity of the flow decreases and materials accretes 
with large speed. Two down-panels of figure |2]show the be- 
havior of sound speed and vertical thickness of the accretion 
disc respectively. As it can be seen from these panels, by 
the increase of magnetic diffusivity (r)o) from — 0.1 the 
isothermal sound speed as well as vertical speed increase. 
Although we see that when the toroidal magnetic field in- 
crease, at first, the diagrams decrease and then increase. It 
means that there is a minimum in two panels. 

In figure [3] we have plotted the coefficients CjS with z- 
component of magnetic field for diffracts values of 770. We 
first find out a strong z-component of magnetic field (/3 Z ) 
leads to an decrease of the infall velocity v r , rotational ve- 
locity isothermal sound speed c s and vertical thickness 
of the disc H/r, which means that a strong magnetic pres- 
sure in the vertical direction prevents the disc matter from 
being accreted, and decreases the effect of gas pressure as 
accretion proceeds. Although it is shown from all four pan- 
els of figure[3] the magnetic diffusivity increases c^s except 
C2 (the rotational velocity). 
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4 Summery and Dissussion 

In this paper we studied the influence of the resistivity on 
the dynamics of advection dominated accretion flows in 
the presence of global magnetic field. Some approxima- 
tions were made in order to simplify the main equations. 
We assumed an axially symmetric, static disc and the a- 
prescription is used for the kinematic coefficient of viscos- 
ity (y = ac s H) and the magnetic diffusivity (77 = ijqc s H). 
Also we ignored the self-gravity and the relativistic effects. 

We have extended Akizuki & Fukue 2006 and Zhang & 
Dai 2008 self-similar solutions to present dynamical struc- 
ture of advection dominated accretion flow by considering 
a three component magnetic field and magnetic diffusivity. 
We have accounted for this possibility by allowing input pa- 
rameters a, rjo, (3 r , f3 v , (3 Z , f to vary in our solutions. 

First, we have shown that the physical quantities of the 
disk are sensitive to the magnetic diffusivity parament i]q. 
We found out the magnetic resistivity increase the radial 
infall velocity v r , the isothermal sound speed c s and the 
vertical thickness of the disk H/r. On the other hand, our 
numerical results shown that when the magnetic diffusivity 
parameter raise up, the rotational velocity of the disk 
decrease. 

Also we have shown that the strong magnetic filed in the 
radial direction increase v r , |i> v | and H/r although c s will 
be decrease which are satisfied with the results were pre- 
sented by Zhang & Dai 2008. Also when the toroidal mag- 
netic field become stronger, our physical quantities v r , \v v \, 
c s and H/r increase (see, e.g., Akizuki & Fukue 2006; Ab- 
bassi et al 2008, 2010; Bu & Yuan & Xie 2009). Moreover 
according to Zhang & Dai 2008, a large vertical magnetic 
field prevents the disk from being accreted and decrease ef- 
fect of the gas pressure. 

Although our preliminary self-similar solutions are too 
simplified, they clearly improve our understanding of the 
physics of ADAFs around a black hole. In order to have 
a realistic picture of an accretion flow a global solution is 
needed rather than the self-similar one. In our future studies 
we intend to investigate the effect of wind and thermal con- 
duction on the observational appearance and properties of a 
hot magnetized flow. 
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Fig. 2 Numerical coefficient as as a function of magnetic parameter f3 v for several values of T]o. The dotted, dashed and solid lines 
correspond to 770 = 0.0, 0.05 and 0.1 respectively. Parameters are set as s — — 0.5(no wind), a — 0.1, f} r = fi z = 1 and f = 1. 
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Fig. 3 Numerical coefficient c^s as a function of magnetic parameter /3 Z for several values of 770. The dotted, dashed and solid lines 
correspond to rjo = 0.0, 0.05 and 0.1 respectively. Parameters are set as s — — 0.5(no wind), a = 0.1, fi r = fi v = 1 and / = 1. 



